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Abstract 

We present additional details of the scanning tunneling microscopy (STM) spectra predicted by the model of pinning 
of dynamic spin collective modes in d-wave superconductor proposed by Polkovnikov et al. (Phys. Rev. B 65, 
220509 (2002)). Along with modulations at the twice the wavevector of the spin collective mode, the local density of 
states (LDOS) displays features linked to the spectrum of the Bogoliubov quasiparticles. The former is expected to 
depend more strongly on an applied magnetic field or the doping concentration. The spin collective mode and the 
quasiparticles are distinct, co-existing, low energy excitations of the d-wave superconductor (strongly coupled only 
in some sectors of the Brillouin zone) , and should not be viewed as mutually exclusive sources of LDOS modulation. 
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1. Introduction 

The observation of a halo of checkerboard modula- 
tion in the local density of states (LDOS) around mag- 
netic field-induced vortex cores by Hoffman et al. [1] 
in optimally doped Bi2Sr2CaCu2 0g+5 has focused at- 
tention on the nature of the spin and charge collective 
modes in the cuprate superconductors. Charge order 
in this region was anticipated theoretically [2,3,4]. By 
'charge order' we mean here periodic spatial modula- 
tions of all local observables invariant under spin rota- 
tions and time-reversal, such as the exchange and pair- 
ing energies, and the LDOS itself. If the modulation is 
primarily on quantities defined on a bond, rather than 
in the charge density on a site, then such order is better 
denoted as 'bond order'; it was this bond order which 
was theoretically anticipated [2,4]. The general connec- 
tion relating such charge/bond order to the 'square' of 
coUinear, non-bipartite, spin correlations [5,6,7,8] (see 
Section 2 below) was used to develop a model [9,10,11] 
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relating the spatial structure and energy spectrum of 
the charge order to the dynamic spin fiuctuation spec- 
trum as measured by neutron scattering [12,13,14,15]. 
There have also been other discussions of competing 
order within vortex cores, including enhanced dynamic 
antiferromagnetism in a spin-gap state [16,17], static 
Neel order [18] , and other models with static spin den- 
sity wave order [19,20,21,22,23,24]. 

More recent measurements of Howald et al. [25] have 
shown that the LDOS modulations in the vortex halo 
also persists in zero magnetic field, albeit with a much 
weaker amplitude: here the pinning of the charge/bond 
order is presumably due to impurities and other imper- 
fections in the Bi2Sr2CaCu2 08+a crystal. Theoretical 
modeling of the STM spectrum has indicated [26,27] 
that the order is likely to be primarily on bond vari- 
ables, as may be derived from a spin-Peierls state [28], 
i.e. a bond order wave [29]. This is in contrast to the 
prediction of a charge density modulation on the Cu 
sites in the anti-phase domain wall picture of 'stripes' 
[6,7,8], which is not observed [25]. 

Hoffman et al. [30] have also reported related mea- 
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surements in zero magnetic field, and have focused at- 
tention on modulations of the LDOS at wavevectors 
throughout the Brillouin zone, and on their disper- 
sion as a function of energy. They have identified some 
wavevectors with quasiparticle excitations of the su- 
perconductor scattered by random impurities, as sug- 
gested in Refs.[31,32]. 

The purpose of the present paper is to discuss the 
relationship between the quasiparticle and collective 
mode contributions to the LDOS modulations in a d- 
wave superconductor. We will use precisely the model 
proposed earlier by Polkovnikov et al. [11] for the pin- 
ning of the spin collective mode: we show that the 
earlier expressions for the induced LDOS modulations 
also contain additional contributions which can be in- 
terpreted as those due to quasiparticle scattering, with 
a spectrum related to observations [30]. The relative 
strength of the collective and quasiparticle contribu- 
tions can be varied by changing the gap to the spin col- 
lective mode: experimentally the gap can be reduced 
(and the collective mode contribution enhanced) by 
lowering the doping or by applying a magnetic field [3], 
and this should allow a test of our predictions. Some re- 
lated observations have also been made independently 
by Han [33]. Our main conclusion is that the LDOS 
modulations observed by STM contain contributions 
from both the spin collective mode and the quasiparti- 
cles: as described in our simple model [11], these are rel- 
atively independent low energy excitations of the BCS 
superconductor, but they can become strongly coupled 
in some sectors of the Brillouin zone if permitted by 
energy and wavevector conservation [34]; they should 
not be viewed as mutually exclusive sources of LDOS 
modulations. 



2. Model 

We assume the cuprate superconductor has a domi- 
nant spin collective mode near the wavevectors K^, = 
(37r/4, tt) and Ky = (tt, 37r/4), and so write for the spin 
operator on the lattice site r at imaginary time r 

5c(r, r) = Re [e"^-'-<I>,„(r, r) + e*'^«--<I-,„(r, r)] , (1) 

where <tx,ya are the spin density wave (SDW) order 
parameters which are assumed to be smooth functions 
of spacetime. As an aside, and for completeness, we 
mention that charge/bond order parameters can also 
be constructed out of the ^x,yoi- We define 

Qa(r,r) s S„(r,r)5a(r a,T) 

« Re [e^^'^-'+^-^^LCr, r)] + . . . (2) 

where a is a vector representing a bond (say the near- 
est neighbor vector), we have assumed that $xa does 



not vary significantly over the spatial distance a, and 
the ellipses denote numerous other similar terms which 
can be deduced from (1). For a a nearest-neighbor vec- 
tor, Qb. measures the modulations in the exchange en- 
ergy, while for a = 0, Qa. measures the local charge 
density (for the t-J model the charge density is lin- 
early related to the on-site S^). Thus we see from (2) 
that, as noted earlier, the charge/bond order param- 
eter is proportional to the square of the SDW order 
parameter. In a similar manner, we will see below that 
the LDOS modulations are also sensitive to the square 
of the SDW order and so occur at wavevectors 2Kx,y: 
however our computation will include retardation ef- 
fects missing from a phenomenological correspondence 
analogous to (2). All our computations below will be 
in a regime where 

(So) = = {^y<.) = (3) 

(because spin rotation invariance is preserved), but we 
will nevertheless have the rotationally invariant 

^ 0, (4) 

a 

because translational invariance is broken by a pinning 
potential. This should be contrasted with the distinct 
physical pictures of Refs. [18,19,20,21,22,23,24] all of 
which break spin rotation invariance with {Sa) ^ 0. 

To keep this paper self-contained, we now briefly de- 
scribe the model used to compute the LDOS modula- 
tions; for details see Ref. [11]. The effective action for 
the $x,ya fluctuations is assumed to be Gaussian: 

5*=T ^ X;r^(q,'^n)|$:ra(q,W«)]^ + (X ^ 3/) (5) 

where ujn are Matsubara frequencies. The Xx.yisi,^) 
denote the dynamic spin susceptibility as a function 
of wavevector q (measured from the ordering wavevec- 
tors 'K.x,y) and frequency uj, which is measured by in- 
elastic neutron scattering. In the absence of available 
experimental input for Bi2Sr2CaCu208-i-5, we will con- 
tinue to use the simple phenomenological form [11] ex- 
pected near a SDW ordering quantum phase transition 
in a superconductor [10]: x~^{'\^'^n) ~ (x)^ + c'^ + , 
where c is a spin wave velocity, and A is the spin gap 
to collective excitations near the wavevectors 'i^x,y (we 
have assumed Xx = Xy = X for simplicity). Han [33] 
has used the same model, but with microscopic com- 
putations on the lattice t-J model as input for x- 

The all-important ingredient which produces an ob- 
servable static modulation in the STM signal is the pin- 
ning of the sliding degree of freedom of the SDW fluc- 
tuations. This pinning is realized by any imperfection 
(impurity, vortex core, etc.) which breaks the lattice 
translational symmetry. We employ the simplest pos- 
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sible pinning term which is invariant under the SU(2) 
spin symmetry, 

.Spin = -Y^jdr [Cx$L(ro, t) + {x ^ y) + ex.] , (6) 

where C^x,y are complex couphrig constants representing 
the pinning potential. 

The electronic quasiparticles of the d-wave super- 
conductor are described by a standard BCS model: 



k 



(7) 



Here = (ckt , cl^^ ) is a Nambu spinor at momentum 
k = {kx, ky), T^'^'^ are Pauli matrices in particle-hole 
space, and /n is a chemical potential. For the kinetic en- 
ergy, Ek we have first (t) and second {t') neighbor hop- 
ping, and the gap function is Ak = (A()/2)(cos fc^, — 
cosA'y). The full Hamiltonian for the conduction elec- 
trons includes a coupling to the collective SDW fluctu- 
ations: 



H = Hbcs + 7/ cl{r)a'i^^c^{r) x 



r 

iK^ -r 



(8) 



+ "I'^«(r)e~ 



+ (a: ^ y)) , 



where a" are the Pauli spin matrices, and 7 describes 
the scattering of the Bogoliubov quasiparticles of the 
superconductor off the SDW fluctuations. 

The predictions of this model [11] for the quasi- 
particle LDOS, to second order in the quasipaxticle- 
coUective mode coupling 7 and to first order in the 
pinning strength (, can be written as 

Sp{r,uj) ^ ImTrii^ X 

Go(r-ri,a;)E(ri,r2,a;)Go(r2 -r,a;), (9) 

where Go is the electronic Green's function in the 
canonical translationally invariant BCS superconduc- 
tor described by (7). There is a contribution to the 
self energy E — E^p from the coupling to the spin 
collective mode, and we focus only on the part which 
breaks translational invariance because of the pinning 
of the collective mode 

SBp(ri,r2,a;) « U^-y^ j ^x(r2, !^)x(ri, i^) 
X Go(r2 - ri,a; - i>) cos(Ka, • (ri -|- r2) - S) 
+ {x^ y). (10) 

The parameter S, which is the phase of the coupling 
C, distinguishes site centered SDWs {d = 0) and bond 
centered SDWs (J = n/2). Following the STM experi- 
ments [1,25,30], we will compute the Fourier transform 
of the LDOS: 



5psp(q,(x') =Tr-!-^^{lm ^ ^^Go(k,w)x 

[Esp (k, k q, tj) 4- Esp(-k, -k - q, w)]Go(k + q, ui) 
cPk 



-i Re 



■Go(k,tj)[Esp(k,k-|-q,u;) 



(2^)2 

-Esp(-k, -k - q, tj)]Go(k + q, a;) } 



(11) 



Note that the imaginary part of this expression van- 
ishes in the presence of inversion symmetry, 5p{r, cj) = 
Sp{—r,u>). The results (10,11) arc for a single pinning 
center, and the realistic many impurity case will be sim- 
ilar provided interference between impurities is small. 

It is not difficult to see from (9-11) that in the limit 
of a large spin gap (A ^ 00), the modulations in the 
LDOS predicted by (11) arc identical to those pre- 
dicted in a model of quasiparticle scattering off iso- 
lated, point-like static impurities or other imperfec- 
tions considered in Ref. [32]. For a general scattering 
potential w(q), and impurity- induced variation in the 
superconducting gap w(qi,q2), the expression for the 
LDOS modulation due to static impurities in the Born 
approximation, (5pimp(q, w), is formally identical to Eq. 
(11) with Esp replaced by 



Sinip(k, k -h q,a;) = u{q)Tz + w(k, k -|- q)Tx . 



(12) 



The following section will compare the predictions of 
(12) to those in (11) of a pinned dynamic SDW with a 
finite A. 

It is interesting to note from (11,12) that for any 
real scattering potential, with 'ft*(q) = «(— q), the 
LDOS modulation at wavevector q is determined by 
J d^fclm [Go(k)Go(k + q,a;)]. Contributions propor- 
tional to J d^klm [Go(k)] Ini [G()(k + q,Ci;)], which are 
used by Ref. [32] to explain the STM observations of 
Refs. [30,35], can only come from lifetime effects which 
produce a complex self energy. 



3. Results 

Sample results from the dynamic SDW model (11) 
are shown in Figs 1-3. For small spin gap A, the LDOS 

shows the expected [5,9] modulations near q — iKx.y 
The magnitude of this modulation depends strongly on 
the value of A, and the strong expected dependence of 
A on doping and applied magnetic field [3] (A decreases 
with decreasing 5 and increasing field) will lead to a 
corresponding variation in the LDOS peaks near q = 
2'Kx,y This peak also disperses as a function of co, with 
the dispersion being stronger for larger values of A as 
shown in Fig. 3; the decrease in the q value of the peak 
with increasing uj is consistent with observations [30] . 

As is clear from Figs. 1 and 2, the LDOS also has 
peaks in its modulation at a number of other wavevec- 



3 




0=30 meV A=10 meV 

Fig. 1. Fourier map of the LDOS 5psp(q, in (11) for A = 10 
meV and quasiparticic dispersion £k 3-nd pairing amplitude 
Ak chosen as in Rcf. [11]. The maxima around (-7r/2,0) come 
mainly from scattering off a pinned SDW. The other peaks are 
due to details of quasiparticle spectrum. 

tors. These arise from the two flanking quasiparticle 
Green's functions in (11), and are consequently sensi- 
tive to all details of quasiparticle spectrum. These ad- 
ditional features are also present in the static impurity 
model (12) whose results are shown in Fig. 4 . It is 
instructive to compare the static impurity scattering 
prediction in Fig 4 to the dynamic spin collective mode 
predictions in Figs. 1 and 2: apart from the strong con- 
tributions at q = i:2K.j:^y present in the latter, the for- 




Fig. 2. As in fig. 1 for A = 40 meV. 



mer has the same general features as the large A cases 
of the latter. This is one of the main points of this pa- 
per. 

The physics of the structure in the static impurity 
scattering LDOS modulation, (5pimp(q, oj) in (12), has 
already been discussed by Wang and Lee [32] . Assum- 
ing spatial inversion symmetry, the expression (12) can 
be rewritten as: 
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Fig. 3. Fourier component of the LDOS, (5p3p(q, 0), for 
q — {qir/a, 0); a large contribution near g — is not shown. 



^Pimp(q,t^) = 2 



jS{-Lo'' +el + Al) X 



(27r)2 

(g; ~ £k)('^ - £k-^q) - AkAk-^q 



k-fq 



+ A 



k+q 



.(k,k + q) ("-"-)f-+'^,+ ^-(V"-+'^) ). (13) 



k+q 



+ a: 



k-l-q 



This expression is not a convolution of the product of 
the density of states of quasiparticles at energy uj with 
momenta differing by q, but rather the product of the 
density of states of one quasiparticle with the real part 
of the Green function of the other. It contains weak 
square- root divergences at special nesting wavevectors 
[36,37], where constant energy contours with momenta 
k and k-l-q [co^ = e^ + Ak = ^k+q + Ak+q) are tangent 
to each other at one of their the points of contact. Such 
wavevectors q form different lines in the Brillouin zone, 
but these are barely visible in our plots. It has been 
claimed [32] that the structure at related wavevectors 
could be enhanced for strong impurity scattering in a 
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Fig. 4. Fourier map of LDOS for scattering off a single localized 
defect, (5pimp {q, tij) in (12). We used a localized defect with 
u(q) — 1, ij{k, k + q) — 2 cos(fca; + Qx) — 2 cos(fcy -t- qy). 

description beyond the Born approximation. Instead, 
in Fig 4, the main contributions for q along the Bril- 
louin zone diagonals come from values of k and k-l-q 
near the nodal points and the underlying Fermi sur- 
face near ±(7r, 0), ±(0, tt) respectively. We emphasize 
that both these features of quasiparticle-scattering in- 
duced LDOS modulations also appear in our model 
[11] of a pinned dynamic spin collective mode. In par- 
ticular, as we noted earlier, if A is large, then the 
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pinned SDW model is equivalent to an effective short- 
range static quasiparticle potential, along with addi- 
tional weak modulations at it^Kx.y (compare Figs 2, 
and 4). For smaller A, the effective quasiparticle poten- 
tial created by the SDW fluctuations is dynamic, but 
the qualitative picture at wavevectors different from 
±2K.x,y remains the same. 

In conclusion, we have shown that the model of 
Ref [11] displays many of the features observed in 
recent STM experiments. Besides the modulations of 
the LDOS at wavevectors ±2Ka;,y arising from the 
spin collective mode, we find other features related to 
the spectrum of the fermionic Bogoliubov quasiparti- 
cles (similar physics is also expected from a possible 
charge/bond order collective mode). The positions of 
the maxima at iiz2K.x,y have a relatively small (but 
not negligible) dependence on the tunnelling energy to, 
while other features are more strongly dependent on 
u and the details of the quasiparticle spectrum . We 
have also made predictions on the doping and mag- 
netic field dependence of the relative contributions 
of the quasiparticle and spin-collective mode induced 
modulations. 

Note added: Some related points and observations ap- 
pear in [38] . 
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